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We present 3D numerical simulation results of moving vortex lattices in presence of ID correlated 
disorder at zero temperature. Our results with field tilting confirm the theoritical predictions of a 
moving Bose glass phase, characterized by transverse pinning and dynamical transverse Meissner 
effect, the moving flux lines being localized along the correlated disorder direction. Beyond a 
critical transverse field, vortex lines exhibit along all their length a "kink" structure resulting from 
an effective static "tin roof pinning potential in the transverse direction. 



Depinning and dynamical ordering of periodic lattices, 
such as vortex lattices (VL) driven by an external force in 
presence of static point disorder has been widely studied 
these last years. A large velocity expansion [l| showed 
that the motion reduces the effects of disorder. In this ap- 
proximation the effects of disorder on moving VL results 
only in an additional "shaking" temperature T s h ~ 1/v 
with expectation of moving crystalline VL at large ve- 
locities. It was however shown in that, due to the 
components of the disorder and displacements transverse 
to the motion, a static disorder persists even at large 
velocity, and leads to a moving glass phase (MG). The 
motion of particles in the MG occurs along rough static 
channels. Depending on the velocity the channels can 
be coupled or decoupled due to dislocations, thus giving 
two generic moving glasses. For weak disorder or large 
velocity in d = 3, channels are elastically coupled leading 
to a topologically ordered structure: the moving Bragg 
glass (MBG) . For larger disorder or in d — 2 channels 
are decoupled by dislocations giving a moving transverse 
glass (MTG) or smectic 0, . The predicted channel 
structure and dynamical ordering was confirmed by nu- 
merical |j, la, Ifl L3l an( i further analytical [8j studies and 
observed experimentally [fj. 

Within the MG model , the situation for correlated 
disorder such as heavy-ion columnar tracks has been re- 
cently investigated 10]. For all velocities at T = or 
for weak velocities at T > it was predicted a moving 
Bose glass phase (MBoG) characterized by a transverse 
critical force and a diverging tilt modulus due to local- 
ization effect arising from the columns. This novel fea- 
ture specific to correlated disorder results in a vanishing 
tilt response below a critical transverse field. This so- 
called dynamical transverse Meissner effect (DTME) is 
a crucial consequence of the MG theory; its observation 
should provide an unambiguous signature of the MBoG 
phase. However, up to now there is no experimental nor 
numerical simulation evidence for this phase. 

This Letter reports 3D numerical simulation of field 
tilting in a moving layered vortex system at zero temper- 



ature in a random pinning landscape of columnar defects. 
The elastically moving vortex lines confirms the existence 
of MBoG displaying DTME. Above a critical angle, the 
vortex lines exhibit a " kink" structure along their length. 
This results from the existence of an effective static pin- 
ning landscape which is almost periodic in the transverse 
direction to the vortex flow. Finite size effects are dis- 
cussed within a simple model of a single flux line in a 
static periodic potential. In particular, it confirms the 
existence of a moving "kink" structure observed in our 
simulation. 

We model a stack of N z Josephson-coupled parallel su- 
perconducting planes of thickness d with interlayer spac- 
ing s. Each layer in the (x, y) plane contains iV,, pancake 
vortices interacting with a random pinning background of 
ID correlated disorder, namely N p columnar pins paral- 
lel to the z direction. At T = the overdamped equation 
of motion of a pancake i in position rj(z) reads 

^ = -E v ^^'^)-E v < cr *Cftp) +F£+FtiIt ( j8 ) 

where pij and Zij are the components of ry = rj — rj 
in cylindrical coordinates, pi P is the in-plane distance 
between the pancake i and a pinning site in the same 
layer at r p , and Vj is the 2D gradient operator act- 
ing in the [x, y) plane. The viscosity coefficient is r\, 
F L = F L x is the Lorentz driving force due to an ap- 
plied current and F tllt (z) is the surface force due to 
the field tilting away from the z axis in the y direc- 
tion. This force acts as a torque on each flux line, i.e. 
F tilt (z = 0) = -F tiH (z = N z s) = F UH y and F ult {z) = 
for pancakes in the bulk. The tilting force modulus is 
defined by F ult = e 2 <foH y / '4tt = 8,Tre 2 e \ 2 ab H y /V3a 2 H z , 
where eo = (</>o/47rA a {,) 2 , \ a b is the in-plane magnetic 
penetration depth, ao is the average vortex distance, 
H y is the transverse field component, and e is the 
anisotropy parameter. The intra-plane vortex-vortex re- 
pulsive interaction is given by a modified Bessel function 
U m ( Pij ) = 2e dK ( Pij /\ ab ). Following 11], the inter- 
plane attractive interaction between pancakes in adjacent 
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layers of altitude z and (z + s) reads U vv (pi 



s) 



(2se /7r)[l + ln(X ab /s)][(pij/2r g y - 1] for Plj < 2r g and 
U m ( Pij , Zij = s) = (2ae /7r)[l + ln{\ ab / s)][ PlJ / r g - 2] 



otherwise; in this expression r g = £ Q b/e, where £ a ;, is the 
in-plane coherence length. Finally, the attractive pinning 
potential is given by U vp {pi P ) = —aA p e~( pi >'> R *'> , where 
A p = (e d/2)ln[l + (i?p/2^ b )] as proposed in 12] and a 
is a tunable parameter. We consider periodic boundary 
conditions of (L x , L y ) sizes in the (x, y) plane while open 
boundaries are taken in the z direction. All details about 
our method for computing the Bessel potential with pe- 
riodic conditions can be found in [13(. Molecular Dy- 
namics simulation is used for N v = 30 vortex lines in a 
rectangular basic cell (L x ,L y ) = (5, 6v / 3/2)A Q b, and for a 
number of layers varying from N z = 19 up to N z = 149. 
The number of columnar pins is set to N p = 30. We 
consider the London limit k = \ a b/£,ab = 90, with an 
average vortex distance ao = \ a b, and d — 2.83 10 _3 A a b, 
s = 8.33 10- 3 A a6 , R p = 0.22 A afc , e = 0.01, 17 = 1. 
These values are coherent with strong anisotropy layered 
superconductors with B z ~ 1000 G. The choice of the 
second order Runge-Kutta algorithm time iteration step 
St is dictated by the dominant force F L , and we take 
St = 7 10~ 2 io, where to = r]^ a b/F L . We choose the tun- 
able pinning parameter a = 1/25 so that the maximum 
pinning force is F^ ax ~ Fq/5 where F = 2e d/X a b is 
the unit force defined by the Bessel interaction. Finally, 
the driving force applied along a principal vortex lattice 
direction x is set to F L ~ 2.8i*o. This corresponds to the 
fully elastic flow limit since F L ~ 25F C L where F C L is the 
critical Lorentz force along x. 

The "experimental" procedure is the folllowing: we 
start in 2D with the (x, y) plane by randomly throw- 
ing N v pancakes and N p gaussian pins, and relaxation 
with zero Lorentz force yields a vortex structure with 
dislocations. The Lorentz force is then slowly increased 
far in the elastic phase up to F L ~ 2.8Fb with steps 
SF L ~ 3.5 10- 3 F every 3.6 10%. The successive dy- 
namical regimes observed in our 2D vortex system are 
the following: (i) pinned regime where all vortices have 
zero velocity; (ii) plastic channels flowing through pinned 
regions and where the motion is periodic (washboard fre- 
quency) as seen in 0, 0] ; (iii) plastic turbulent flow with 
no stationnary vortices as seen in [6J, |7j ; (iv) weakly de- 
coupled channels named MTG since quasi long range lon- 
gitudinal order exists in each channel; (v) elastic phase 
where all vortices have the same average velocity and 
named MBG since no dislocation appears and the mo- 
tion occurs through rough static channels even in the high 
driving phase F L ~ 25F C L studied below. Our 2D numer- 
ical results are therefore in agreement with the moving 
glass theory developed in Q • Finally, note that they are 
consistent with the phase diagram region corresponding 
to the intermediate pinning of Ref 0. 

Now starting from this situation, we extend the vor- 
tices and the pins in the z direction in order to obtain 3D 
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Figure 1: Average vortex line inclination tan 9b = B y /B z 
versus the field inclination tandn = H y /H z for several sys- 
tem sizes N z = 19,29,39,49,149 (filled circles). The linear 
response of the pinning free vortex system is shown with open 
circles. Size effects show that the partial transverse screening 
observed for H y /H z -< 0.50 tends to be total in the infinite 
size limit N z — » 00, showing the existence of the DTME below 
the critical transverse field Hy/H z ~ 0.50. Insert: projection 
in the (y, z) plane of the 30 vortex lines composed of N z = 149 
layers illustrating the MBoG phase (thick lines) and DTME, 
compared to the vortex lines of the non-disordered system 
(thin oblique lines). In the MBoG phase, the transverse field 
penetrates only at the surface and vortex lines in the bulk are 
aligned with the columnar defect direction z. In the pinning 
free situation, the vortex lines form straight lines in an oblique 
direction defined by the field inclination and Josephson cou- 
pling. Note that the whole vortex lattice is moving elastically 
in the x direction, i.e. perpendicularly to the image. 



straight vortex lines moving elastically along the x direc- 
tion through a random columnar pin landscape. Then, 
starting from this situation, we now slowly tilt the mag- 
netic field H away from z in the y direction by increasing 
the H y component transversally to the vortex flow. For 
several system sizes, Fig. ^ displays the vortex line re- 
sponse to low field tilting, i.e. the average vortex line 
inclination tan^s = B y /B z versus the field inclination 
tan(9# = H y /H z , compared with the linear response of 
the pinning free vortex system. It shows a partial screen- 
ing of the transverse component H y of the field, which 
tends to be total for infinite vortex line length since the 
slopes tend towards zero when N z — * 00. This trans- 
verse screening is followed by a transition at a critical 
field H y /H z ~ 0.50, above which the average vortex line 
inclination recovers the one obtained without any pin- 
ning. Therefore, Fig. ^ shows the existence of DTME 
which is the very signature of the MBoG phase as pre- 
dicted in 01 • a good illustration of DTME, we show 
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Figure 2: Plateaux and hysteresis observed above the transi- 
tion for N z —49. The origin of each plateau is connected to 
a given number of "pinned" regions passed through by each 
vortex line. For details about the plateaux labels and the 
"pinned" regions defining the effective static pinning land- 
scape, see Fig. Eland text below. Finite size effects show that 
the plateaux and hysteresis disappear in the infinite size limit 
N z — > oo. Plateaux and hysteresis may therefore be observed 
in thin real samples. 

in the insert of Fig. ^ the vortex line configuration of 
the MBoG compared with the non-disordered lattice ob- 
tained for the same field tilting. One clearly sees that 
below the transition the transverse field enters at the 
surface without penetrating into the bulk, in opposition 
to the pinning free vortex system. So below the critical 
transverse field, the vortex lines remain parallel to the 
columnar defects even though they are moving at enough 
large velocity. 

We shall now examine the situation above the crit- 
ical transverse field. Far above Hy/H z and whatever 
the vortex line size is, the average vortex line inclina- 
tion recovers the one obtained without any pinning (not 
shown). A more interesting case happens in the inter- 
mediate transverse field range above the transition and 
shown for N z = 49 layers in Fig. El Increasing the trans- 
verse field above its critical value Hy/H z ~ 0.50 gives 
rise to plateaux in the tan#s vs H y /H z curve. A con- 
sequence of these plateaux is that decreasing the trans- 
verse field generates hysteresis in the vortex inclination. 
A typical vortex line configuration above the transition 
is shown for N z — 149 layers in Fig. |3fa). A modulated 
structure z(y) appears along their length. This modu- 
lation is more obvious looking at the derivative dz/dy 
shown for all the vortex lines in Fig. E^b) . By using the 
minima and maxima of the derivative we may roughly 
define two parts of small and big inclinations in the vor- 




Figure 3: (a) Projection in the (y, z) plane of the 30 mov- 
ing vortex line composed of N z = 149 layers and obtained 
for Hy/Hz ~ 0.53 while increasing the transverse field. In 
grey are materialised the "pinned" regions which define an 
effective pinning landscape which is static and almost peri- 
odic, (b) Derivative dz/dy for all the 30 vortex lines shown in 
(a). The vortex line modulations along their length become 
clearly visible. Taking into account the minima and inflex- 
ion points of the derivative allows the determination of the 
"pinned" regions. 

tex lines with respect to z. The small ones are termed 
"pinned" regions and are linked by "unpinned" regions 
where the vortex inclination is bigger. The materiali- 
sation of the "pinned" regions is shown in grey in Fig. 
|3Ja) . We visualize in this way a static pinning landscape 
in the transverse direction y with effective columnar pin 
width corresponding to twice the input pinning "radius" 
R p . Each plateau of Fig. |2Jcan be thereby labelled with 
an odd number corresponding to the number of " pinned" 
regions passed through by each vortex line. In this way, 
while increasing the transverse field above the transition, 
the vortex lines firstly pass through 3 pinned regions for 
a given range of H y defining the "P3" plateau, then the 
system jumps to the following plateau "P5" where the 
vortex lines pass through 5 pinned regions, and so on. 
The same process holds while decreasing the transverse 
field generating thereby hysteresis. Finally, the central 
point displayed in Fig. 01a) is that the vortex lines which 
are moving in the x direction perpendicularly to the fig- 
ure plane see an effective static tin roof potential in the 
transverse direction y. This is a direct consequence of 
the existence of static pinned channels of motion 0, ^| , 
and clearly shows that the action of the disorder in the 
moving system cannot be simply reduced to a tem- 
perature. Note that this effective potential is not strictly 
periodic because of the channel roughness. 
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We now turn to the discussion of the finite size effects 
in our simulation by considering the tilt response of a 
single static elastic vortex line in a tin roof pinning po- 
tential at T — 0. Such a simple model allows to get 
a valuable insight into the transverse properties of the 
moving vortex lattice in presence of correlated disorder, 
e.g. the TME and the kink vortex line structure. The 
energy of a line of length L is given by: 

E{v) = f g dz {l{^) + V e f f (u)\+f(u(L)-u(0)),( 

where u(z) is the one-dimensional displacement field in 
the y direction, c = e 2 eo is the elastic constant, / oc H y 
is a surface force, and V e ff(u) = g [l — cos^Tra/VSao)] 
with g being a constant. There exists a ground state ulj 
which minimizes E(u). We calculate the tilt response of 
the vortex line defined by tan 9 = (ulj{L) — ulj{0))/L. 
In the limit L — > oo, we find a critical force f c — A^fgc/ir 
below which u_l =00 ./ = 0, giving tan (9 = and hence 
TME. For f > f c the vortex line is tilted and its ground 
state ul=ooj{z) are expressed in terms of Jacobian ellip- 
tic functions which generate kink vortex line structures 
as observed in our simulation (see Fig.^a)). The tilt re- 
sponse is shown in the insert of Fig. 0] together with the 
pure elastic response, i.e. g = 0. In the latter case, the 
ground states are linear functions and hence the vortex is 
a straight line of slope tan 9 — f /c. The kink vortex line 
structure obviously rely only on the periodic potential. 
Note that /* = n/3//8 7rc is the equivalent quantity of 
H y /H z . For a finite size line and / — 0, the minimum 
energy solution is ulj=o = 0, i.e. the line is pinned. If / 
is nonzero but sufficiently weak, both extremities of the 
line are pulled away from its minimum energy configura- 
tion in a continuous way, yielding a nonzero tilt response 
as observed in our simulation (see Fig.QJand insert). Fur- 
thermore, for given values of /, the vortex line abruptly 
stretches to a more favourable position resulting in jumps 
in the tilt response. We show in Fig. 0] the tilt response 
computed for L = 0.4 ao which is the equivalent size of 
the 49 layers in our simulation results shown in Fig. |2 
the straight line displays the pure elatic response g = 0. 
We therefore find plateaux and hysteresis which are very 
close to the ones observed in our simulation and shown 
in Fig. □ 

In conclusion, our numerical simulation results give the 
first confirmation of the existence of the MBoG phase dis- 
playing a DTME developed in [lo|. Furthermore, above 
the critical transverse field, the simulation revealed a 
"kinked" vortex line structure resulting from an almost 
periodic effective static pinning potential in the trans- 
verse direction. 
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Figure 4: Plateaux and hysteresis observed within the simple 
model (Eq. for a finite length L = 0.4 ao and g/c= 18.75. 
Black filled circles represent the metastable states followed by 
the elastic line when the surface force / is slowly increased, 
then decreased. The white filled circles represent the true 
equilibrium states. Insert: Eq. Q results in the infinite size 
limit L — > oo showing the existence of a critical dimension- 
less force /c = \/3/27r 2 .\/g/c and the disappearance of the 
plateaux and hysteresis. 
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